Turbulence Modeling Approaches Research Articles

Wall-modeled large eddy simulation of a tandem wing configuration in transonic flow

In this study, the interaction between a turbulent wake and the boundary layer of a horizontal tail plane (HTP) in the transonic flow regime is investigated. The setup considered corresponds to a generic tandem wing configuration with an OAT15A airfoil as the main wing and a National Advisory Committee for Aeronautics (NACA) 64A-110 as an HTP. Due to the transonic flow, the suction side of the OAT15A exhibits buffet. The numerical approach consists of a two-stage procedure in which a detached eddy simulation (DES) provides unsteady inflow conditions for a subsequent zonal high-fidelity large eddy simulation (LES) performed for the HTP region only; the turbulent boundary layer is modeled using scale-resolved wall-modeled LES (WMLES). The study mainly pursues two objectives: first, to discuss the influence of wake turbulence on the flow characteristics of the NACA airfoil; and second, to evaluate the capability of WMLES as a low-cost but high-resolution numerical approach in challenging flow conditions. Essentially, the study confirms the expected result that the outer part of the HTP's boundary layer is dominated by the wake of the main wing. Mainly based on a discussion of turbulence spectra, the study further demonstrates the advantage of WMLES over DES, proving that WMLES is able to capture the effects of wake turbulence on boundary layer dynamics, and thus validates the WMLES approach as a cost-effective, high-resolution turbulence modeling approach. On a superordinate level, the study further sketches a possible way on how a flow problem with such a strongly unsteady behavior could be systematically evaluated.

Improved self-adaptive turbulence eddy simulation for complex flows and stall prediction using high-order schemes

It is challenging to apply numerical simulations to accurately predict the stall behavior of aircraft equipped with high-lift devices. Simulations with Reynolds-Averaged NavierStokes (RANS) models suffer from lack of the reliability at high angles of attack with separated and reattached boundary layers, whereas wall-resolved Large Eddy Simulations (LES) of wall-bounded flows at high Reynolds numbers currently costs too much computational resources. A new unified hybrid turbulence modeling approach, denoted the Self-Adaptive Turbulence Eddy Simulation (SATES), is proposed and applied for complex turbulent flows combining with high-order numerical scheme of the Weighted Compact Nonlinear Scheme (WCNS) in the present study. It enables a seamless evolution from unsteady RANS to LES and finally approaches Direct Numerical Simulation (DNS) depending on the turbulent scales. In the framework of SATES, a new SATES-σ model with an adaptive model coefficient is developed by extending the underlying LES mode based on an enhanced sub-grid-scale model of the σ-model. The new SATES-σ is first examined in two benchmark cases of channel flow and flow past a square cylinder. Then, it is validated in supercritical flow past a circular cylinder to assess the performance of turbulent models. The results show significant improvements over the previous SATES and IDDES in the predictions of boundary layer flow. Finally, successful application is achieved in the accurate prediction of the stall of the MD-30P30N airfoil at a Reynolds number of 9×106 with wide angles of attack. The simulation results show good agreement with experimental results for surface pressure even for the challenging cases of 21 and 23 deg angles of attack. Again, the SATES-σ shows better results than the previous SATES and IDDES. The presented method has considerable potential for the challenging stall predictions.

Modeling of the influence of local heat sources on a light gas stratification formation and erosion in a large-scale experimental facility using eddy resolving numerical approach

The main objective of this work is to evaluate the results of application of Implicit Large Eddy Simulation (ILES) turbulence modeling approach to simulate the complex hydrogen safety problem. The paper presents the results of numerical simulation of the OECD/NEA HYMERES (HYdrogen Mitigation Experiments for Reactor Safety) HP2 series tests in the PANDA facility, aimed at studying the thermal effect of passive autocatalytic recombiner (PAR) operation on mixing the containment atmosphere, using CABARET-SC1 code. The code is based on the eddy resolving CABARET technique, which allows implicit modeling of the subgrid turbulence scales without using tuning parameters (ILES approximation). The natural convection flow inside the PAR model was calculated explicitly using a simplified model of the heating section and the housing. Good agreement of calculated transient with experimental data was obtained, which indicates the effectiveness of such approach.

Towards the understanding of vortex breakdown for improved RANS turbulence modeling

A triple-delta wing configuration in transonic conditions has been investigated numerically, explicitly focusing on the leading-edge vortices during sideslip flow. Several shock waves manifest over the fuselage, leading to vortex-shock interactions. The primary focus of the analysis centers on the vortex breakdown phenomenon and its underlying causes. The flow field obtained from different turbulence modeling approaches is compared. This study aims to identify the limitations of a standard one-equation model in dealing with such complex applications. It explores physical and modeling reasons leading to the model's inaccuracies. The impact of turbulence treatment on the numerical results is explored, discussing turbulence-related variables. A novel and straightforward modification to the one-equation turbulence model is then introduced and assessed compared to experimental and computational data. The proposed model appears promising for enhancing the accuracy of the one-equation Reynolds Averaged Navier-Stokes outcomes.

SIMULATING SPRAY DYNAMICS WITH A FINITE ELEMENT METHOD FOR INTERNAL COMBUSTION ENGINES USING LARGE EDDY SIMULATIONS

Spray dynamics in an internal combustion engine is comprised of complex phenomena while interacting with unsteady turbulence. Physics requires detailed modeling of the dynamics for spray and carrier gases to accurately predict a spray's fate. Large eddy simulation (LES) turbulence modeling approaches are capable predictors of the turbulent processes and of dynamically modeling sub-grid scales, therefore enabling calculation of model coefficients for the smallest resolved scale. Dynamic LES methods are also well suited for unsteady flows associated with spray injection and engine fluid dynamics. In this study, for the first time, a dynamic Verman LES scheme employed in a stabilized a finite element (FE) framework is used to model the spray process with emphasis on injected fuels for simulating internal combustion engines. Spray modeling often comprises a coupled Eulerian-Lagrangian approach to capture the droplet/particle dynamics, where the droplets are modeled in the Lagrangian frame, interacting with the gas. The momentum and heat exchange between the fluid gases and the evaporating and atomizing spray droplets are modeled in a two-way coupling system as described in this paper. Direct injected liquid is modeled in this paper as a spherical ligament of fuel and ligament breakup to atomization use our version of the Kelvin Helmholtz break-up scheme. Discussed are models and methods of the whole system in some detail. The method for simulation of the fluid's momentum, heat transfer, and turbulence are discussed as is the system to evaluate droplet or ligament properties. Validation or results of the modeling are presented on test cases as determined by the Engine Combustion Network, base cases G along with some benchmarking cases. We demonstrate in the paper by using the FE method combined with the Verman LES model, and our version of the Kelvin Helmholtz-Rayleigh Taylor spray break-up model combined with two-way coupling is quite capable of modeling sprays, having distinct and significant advantages over the more traditional finite volume discretization framework using similar spray modeling methods.

Numerical Investigations of Turbulent Flow Through A 90-Degree Pipe Bend And Honeycomb Straightener

Abstract Pipe bends are commonly used in piping systems in offshore and subsea installations. The present study explores the design considerations for the honeycomb straightener inserted downstream of a 90-degree pipe bend. The objective of the study is to evaluate the effectiveness of the honeycomb in suppressing the flow swirling for different distances from the bend outlet (Lb) and different values of the honeycomb thickness (t). The turbulent flow through the 90-degree circular pipe bend with the honeycomb straightener is investigated by carrying out numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) turbulence modeling approach. The Explicit Algebraic Reynolds Stress Model (EARSM) is adopted to resolve the Reynolds stresses. The honeycomb thickness to pipe diameter ratio (t/D) is varied between 0.1 and 1. The normalized distance from the bend outlet to the honeycomb straightener (Lb/D) is varied between 1 and 5. The disturbance in the velocity field is generated by the pipe bend with the curvature radius to pipe diameter ratio (Rc/D) of 2 and Reynolds number (Re) of 2×105. It is found that both the increase in Lb/D and t/D improves the performance of the device in removing the swirl behind the bend outlet. The best performance is observed for the honeycomb straightener with the distance Lb/D=5 and thickness t/D=0.5.

Evaluation of the overset method applied to an active lift turbine

In the current context of energy transition, the development and improvement of energy recovery systems is essential. These system are numerous and some of them use turbines to convert flow kinetic energy into electricity. For a long time, vertical-axis turbines were considered less efficient than horizontal-axis ones. However, they are still better suited to urban areas. Thus, they are now being studied again and new turbine concepts are emerging that can works differently from today’s turbines. The entire dynamics of these new-generation turbines can be completely different from those of conventional turbines. The active lift turbine, for example, changes the trajectory of the blades (by varying the turbine radius) to convert the normal forces applying on them into torque. This particularity hinders the use of conventional rotation methods for turbines (Arbitrary Mesh Interface (AMI), Rotating Frame of Reference (RFR), …) in Computational Fluid Dynamics (CFD) and an innovative method such as overset could be an alternative. The purpose of this work is to evaluate the overset method for these new-generation turbines modelling. It enables the motion of a structure in its environment to be controlled at will. Various modelling parameters are being tested. The impact of the mesh size (mesh convergence), the turbulence modelling approach (RANS/LES) and the interpolation method between the background mesh and the overset mesh (Cell Volume Weight/Inverse distances) are evaluated and compared. The overset method is first validated on a conventional turbine (Darrieus turbine). The results obtained for the conventional turbine are close to those expected, with a relative error (numerical/experimental) between 3% and 8% depending of the Tip Speed Ratio (TSR) with an averaged error around 7%. Then overset method is extended to an active lift turbine to assess the potential of this new technology. The overset method proves its worth for this type of problem, which cannot be solved by conventional CFD methods of setting structures into motion. It provides control over blade motion without loss of information. Overall, this study demonstrates the interest and applicability of the overset method to the Active Lift Turbine or other kind of turbines that have a specific motion.

Data-driven turbulence modelling of inherently unsteady flow in stratified water storage tanks

The flow and thermal characteristics of thermal stratification appears in a wide range of industrial applications related to energy storage. In a previous study (Xu et al. submitted to Flow, Turbulence and Combustion), it has been shown that linear closure models in unsteady Reynolds averaged Navier–Stokes (RANS) calculations are unable to accurately capture the transient flow and thermal characteristics occurring in stratified water storage tank configurations. In the current study, we extend a data-driven framework to develop turbulence and heat flux closures for a three-dimensional (3D) geometry of a stratified water storage tank, using data from a highly resolved large eddy simulation. The transient high-fidelity data set is partitioned into several time-intervals used for closure development and testing. In detail, we first present the model formulations for 3D buoyant flow. We then use the reference data to extract the most suitable input features, i.e. variables used to regress the data driven model. In order to reduce the complexity of the resulting models and accelerate the training process, a recursive basis function elimination is performed, using alignment of the different basis functions with the target turbulence stress and heat flux as metric. Using the identified features and the subsets of basis functions, novel closures are then developed using a symbolic regression framework, based on Gene Expression Programming (GEP). Models are developed either for each time interval or for all intervals at once. The different machine-learnt turbulence and heat-flux closure models are then tested for two different tank configurations where the thermal stratification is affected by positively/negatively buoyant jets within simple or complex three-dimensional geometries. It is demonstrated that the extended GEP-based data-driven turbulence modelling approach is able to create models that can accurately predict flow with inherent unsteadiness in three dimensional complex-geometry cases.

Secondary flow and thermal field dynamics in a T-junction with an upstream elbow: A POD-aided, eddy-sensitized Reynolds-stress modeling study

The thermal mixing of the flow-crossing streams in a T-shaped pipe junction configuration with the main branch inflow forming an elbow is presently simulated by applying a RANS-based eddy-resolving model of turbulence. The low-temperature main stream enters the T-junction with a Reynolds number of Rem=107893 after passing through a 90o-turned elbow located upstream. The high temperature branch stream is injected perpendicularly into the main stream at a Reynolds number of Rem=5422. A significant difference between the flow velocities of the two streams leads to the formation of a characteristic mixing process resembling a reattached jet configuration. Of particular interest is the interference of the branch flow with the main flow characterized by secondary vortices – the so – called Dean vortices, that are convected from the curved main inlet generating a low-rank mixing structures in its wake. The reference LES-based (Large-Eddy Simulation) database is provided by Tunstall et al. (2016). The computational framework relying on the scale-adaptive turbulence modeling approach, introduced principally by Menter and Egorov (2010), is used to simulate the dynamics of large-scale turbulent structures with moderate spatial and temporal resolution requirements. The motion of unresolved subscale structures is described by an appropriately sensitized differential near-wall Reynolds stress model according to Jakirlić and Maduta (2015). In addition, Proper-Orthogonal Decomposition (POD) is used to analyze the prevailing flow and thermal field mechanisms. The results obtained show a reasonable level of agreement with the LES reference, both from the standpoint of first-order statistics and spectral reproduction of the dynamics of coherent structures.

Development of an effective two-equation turbulence modeling approach for simulating aerosol deposition across a range of turbulence levels

Pharmaceutical aerosol systems present a significant challenge to computational fluid dynamics (CFD) modeling based on the need to capture multiple levels of turbulence, frequent transition between laminar and turbulent flows, anisotropic turbulent particle dispersion, and near-wall particle transport phenomena often within geometrically complex systems over multiple time scales. Two-equation turbulence models, such as the k−ω family of approximations, offer a computationally efficient solution approach, but are known to require the use of near-wall (NW) corrections and eddy interaction model (EIM) modifications for accurate predictions of aerosol deposition. The objective of this study was to develop an efficient and effective two-equation turbulence modeling approach that enables accurate predictions of pharmaceutical aerosol deposition across a range of turbulence levels. Key systems considered were the traditional aerosol deposition benchmark cases of a 90-degree bend (Re=6,000) and a vertical straight section of pipe (Re=10,000), as well as a highly complex case of direct-to-infant (D2I) nose-to-lung pharmaceutical aerosol delivery from an air-jet dry powder inhaler (DPI) including a patient interface and infant nasal geometry through mid-trachea (500<Re<7,000). Of the k−ω family of models, the low Reynolds number (LRN) shear stress transport (SST) approach was determined to provide the best agreement with experimental aerosol deposition data in the D2I system, based on an improved simulation of turbulent jet flow that frequently occurs in DPIs. Considering NW corrections, a new correlation was developed to quantitatively predict best regional values of the y+limit, within which anisotropic NW turbulence is approximated. Considering EIM modifications, a previously described drift correction approach was implemented in pharmaceutical aerosol simulations for the first time. Considering all model corrections and modifications applied to the D2I system, regional relative errors in deposition fractions between CFD predictions and new experimental data were improved from 19–207% (no modifications) to 2–15% (all modifications) with a notable decrease in computational time (up to ∼15%). In conclusion, the highly efficient two-equation k−ω models with physically realistic corrections and modifications provided a viable, efficient and accurate approach to simulate the transport and deposition of pharmaceutical aerosols in complex airway systems that include laminar, turbulent and transitional flows.

A Generic Model for Benchmark Aerodynamic Analysis of Fifth-Generation High-Performance Aircraft

This paper introduces a generic model for the study of aerodynamic behaviour relevant to fifth-generation high-performance aircraft. The model design is presented, outlining simplifications made to retain the key features of modern high-performance vehicles while ensuring a manufacturable geometry. Subsonic wind tunnel tests were performed with force and moment balance measurements used to develop a database of experimental validation data for the platform at a freestream velocity of 20 m/s. Numerical simulations are also presented and validated by the experiments and further employed to ensure the vortex behaviour is consistent with contemporary high-performance platforms. A sensitivity study of the computational predictions from the turbulence modelling approach is also presented. This geometry is the first in a suite of representative aircraft geometries (the Sydney Standard Aerodynamic Models), in which all geometries, computational models, and experimental data are made openly available to the research community (accessible via this link: https://zenodo.org/communities/ssam_gen5/) to serve as validation test cases and promote best practices in aerodynamic modelling.

Self-Adaptive Turbulence Eddy Simulation with a high-order finite differencing method for high Reynolds number complex flows

Accurate prediction of high Reynolds number wall-bounded flows is a challenge for turbulence modeling and numerical discretization. Moreover, complex engineering applications increase the difficulties in implementing robust and low-dissipation numerical schemes. A newly developed unified turbulence approach, Self-Adaptive Turbulence Eddy Simulation (SATES), is assessed in the simulation of a multi-element 30P30N airfoil at a high Reynolds number of 1.71 × 106. To reduce the numerical error, high-order finite differencing-based weighted compact nonlinear schemes (WCNSs) are successfully applied. The effects of grid resolution, numerical schemes and turbulence model are investigated by grouping studies. The SATES shows low sensitivity to grid resolution and predicts satisfactory results even on a coarse grid containing approximately 7.2 million cells. Compared with second-order accurate central differencing, the high-order accurate WCNS has less numerical dissipation and significantly improves the resolution of turbulence structures and the accuracy of numerical results. Z-type nonlinear weighted averaging interpolation, designed for shock capturing property in the WCNS scheme, leads to extra numerical dissipation and affects the accuracy of the results. The numerical results of the SATES are more accurate than those of the Smagorinsky Large Eddy Simulation (LES) and the Improved Delayed Detached Eddy Simulation (IDDES). The SATES provides sufficient modeling of the boundary layer on the coarse grid and a fast ‘transition’ from the RANS to the LES mode for accurate prediction of the separated shear layer and flap boundary layer separation. These results confirm that the SATES turbulence modeling approach is efficient for high Reynolds number flows and that the high-order accurate WCNS schemes show sufficient robustness and accuracy for engineering applications.

Mesh-Free Analysis of a Vertical Axis Wind Turbine Using Lattice Boltzmann Method and Various Turbulence Models

This study aims to investigate the aerodynamic analysis of a Darrieus-type vertical axis wind turbine (VAWT) using the Lattice Boltzmann Method (LBM). The objective is to assess the accuracy and performance of the meshless LBM approach in predicting torque coefficients, velocity, turbulence intensity, and vorticity distributions for VAWT aerodynamic analysis. Two turbulence modelling approaches, Large Eddy Simulation (LES) and Reynolds-Averaged Navier-Stokes (RANS), are employed to model the flow domain. The central problem revolves around comparing the performance of different turbulence models based on their agreement with experimental results for power and torque coefficients. The findings demonstrate the effectiveness of the WALE turbulence model in achieving the best agreement with experimental data. Overall, the study provides valuable insights into applying LBM in VAWT aerodynamic analysis and highlights the advantages of the meshless approach compared to traditional CFD methods.

RANS-DES Comparison on the Numerical Modelling of a Downforce Generating Wing Moving over a Solid Surface

Present paper covers the numerical modelling of a 3D wing operating near a solid surface using Reynolds Averaged Navier Stokes (RANS) and Detached Eddy Simulation (DES). The characteristics of the flow over a wing near a surface is significantly different from the free flight. Capturing the changing aerodynamics of an object moving close to a surface have importance for flight and vehicle aerodynamics. This study aims to investigate the reliability of RANS and DES turbulence modelling approaches on the numerical solution of the challenging ground effect flow physics of the downforce generating (or inverted) wings. For this purpose, a 3D symmetrical wing moving in the close proximity of a solid surface is numerically modelled in fully turbulent flow conditions. The results show that, DES approach predicts the lift force slightly better than RANS. On the other hand, tip leakage flow capturing and the separating vortex field modelling of the DES approach is clearly superior. Comparison of velocity, vorticity and the vortex visualization results are presented for the better understanding of the performance of the two modelling approaches.

Direct Numerical Simulation of Low and Unitary Prandtl Number Fluids in Reactor Downcomer Geometry

Mixed convection of low and unitary Prandtl fluids in a vertical passage is fundamental to passive heat removal in liquid metal and gas-cooled advanced reactor designs. Capturing the influence of buoyancy in flow and heat transfer in engineering analysis is hence a cornerstone to the safety of the next-generation reactor. However, accurate prediction of the mixed convection phenomenon has eluded current turbulence and heat transfer modeling approaches, yet further development and validation of modeling methods is limited by a scarcity of high-fidelity data pertaining to reactor heat transfer. In this work, a series of direct numerical simulations was conducted to investigate the influence of buoyancy on descending flow of liquid sodium, lead, and unitary Prandtl fluid in a differentially heated channel that represents the reactor downcomer region. From time-averaged statistics, flow-opposing/aiding buoyant plumes near the heated/cooled wall distort the mean velocity distribution, which gives rise to promotion/suppression of turbulence intensity and modification of turbulent shear stress and heat flux distribution. Frequency analysis of time series also suggests the existence of large-scale convective and thermal structures rising from the heated wall. As a general trend, fluids of lower Prandtl number were found to be more susceptible to the buoyancy effect due to stronger differential buoyancy across the channel. On the other hand, the effectiveness of convective heat transfer of the three studied fluids showed a distinct trend against the influence of buoyancy. Physical reasoning on observation of the Nusselt number trend is also discussed.

Learning Neural-Network-Based Turbulence Models for External Transonic Flows Using Ensemble Kalman Method

This paper presents a neural-network-based turbulence modeling approach for transonic flows based on the ensemble Kalman method. The approach adopts a tensor-basis neural network for the Reynolds-stress representation, with modified inputs to consider fluid compressibility. The normalization of input features is also investigated to avoid feature collapsing in the presence of shock waves. Moreover, the turbulent heat flux is accordingly estimated with the neural-network-based turbulence model based on the gradient diffusion hypothesis. The ensemble Kalman method is used to train the neural network with the experimental data in velocity and wall pressure due to its derivative-free nature. The proposed framework is tested in two canonical configurations, that is, two-dimensional transonic flows over the RAE2822 airfoils and three-dimensional transonic flows over the ONERA M6 wings. Numerical results demonstrate the capability of the proposed method in learning accurate turbulence models for external transonic flows.

Numerical simulation of bluff body turbulent flows using hybrid RANS/LES turbulence models

Many engineering applications involve turbulent flows around bluff bodies. Because of their intrinsically unsteady dynamics, bluff body characteristic flows feature unique turbulence-related phenomena, which makes their numerical modeling challenging. Accordingly, accounting for a circular bluff body flow configuration, three different turbulence modeling approaches are investigated in this work, (i) Reynolds-averaged Navier–Stokes (RANS), (ii) large eddy simulation (LES), and (iii) hybrid RANS/LES. Regarding the hybrid approaches, two variants of the detached eddy simulation (DES) one, delayed DES (DDES) and improved delayed DES (IDDES), are studied. As RANS model, the $$\mathrm{k}-\mathrm{\omega SST}$$ is utilized here. This RANS model is also used as the background one for both DDES and IDDES. Wall-adaptive local eddy viscosity (WALE) is used in turn as the sub-grid scale (SGS) model for LES. The velocity two-point correlation function is used to assess the mesh size requirements. When compared to experimental data, the obtained numerical results indicate that RANS overestimates the recirculating bubble length by over 18% and is not capable of describing the turbulent kinetic energy and the flow anisotropy in agreement with the experimental data. In contrast, LES, DDES, and IDDES are all within 1% of the recirculating bubble length while predicting both the Reynolds stress tensor components and the corresponding flow anisotropy in agreement with the measurements. Besides, normalized anisotropy tensor invariants maxima in the shear layer were reproduced by all scale resolving models studied here, but they failed to yield the local extrema measured within the wake recirculation region. A comparative analysis of the anisotropic Reynolds stress tensor invariances underscores the adequacy of the scale resolving models.

Multiscale model for blood flow after a bileaflet artificial aortic valve implantation

Cardiovascular diseases are the leading cause of mortality in the world, mainly due to atherosclerosis and its consequences. The article presents the numerical model of the blood flow through artificial aortic valve. The overset mesh approach was applied to simulate the valve leaflets motion and to realize the moving mesh, in the aortic arch and the main branches of cardiovascular system. To capture the cardiac system’s response and the effect of vessel compliance on the outlet pressure, the lumped parameter model has been also included within the solution procedure. Three different turbulence modeling approaches were used and compared — the laminar, k-ϵ and k-ω model. The simulation results were also compared with the model excluding the moving valve geometry and the importance of the lumped parameter model for the outlet boundary condition was analyzed. Proposed numerical model and protocol was found as suitable for performing the virtual operations on the real patient vasculature geometry. The time-efficient turbulence model and overall solving procedure allows to support the clinicians in making decisions about the patient treatment and to predict the results of the future surgery.

Zonal Turbulence Modeling Approach for Simulating Compartment Fire for Initial Hazard Assessment

The flow field driven by a compartment fire usually contains several flow zones with different physical structures. As each type of turbulence model has its own predominant application area, it is logical to apply two or more simple turbulence models to the same fire-induced flow field at different locations according to their predominant features to yield a comparatively simple, accurate, and stable zonal turbulence model. A zonal turbulence model, which is a hybrid of the standard k-ε model and its modification, is developed in this paper. The model is tested and compared with the experimental data. A promising improvement is observed when comparing it with the base turbulence model, i.e., the standard k-ε model, especially in the recirculating region near the corners of the compartment. This approach in having different zones in the plume region will be useful for handling more scenarios at the initial stage of fire hazard assessments.

Numerical simulation of scour at the rear side of a coastal revetment

This paper presents the results of a numerical modeling study on the scouring of unprotected rear side material of a rubble mound coastal revetment due to the overtopping of solitary-like waves utilizing a coupled hydro-morphodynamic computational fluid dynamics (CFD) model. Three cases having various wave heights are tested with six different turbulence models together with different wall functions. The hydrodynamic results (free-surface elevations, overtopping volumes, and jet thicknesses) and morphological results (scour profile and scour hole properties) are compared with the existing experimental data to evaluate the performances of different turbulence models and wall functions. It is found that the numerical model can predict the hydrodynamic and morphological results with a high level of accuracy when the appropriate turbulence model and wall function are used. The importance of turbulence in the scouring processes and the turbulence modeling approach in hydro-morphodynamic CFD simulations are highlighted.